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For understanding an anomalous nuclear effect experimentally observed for the beryllium-9 nu- 
cleus at the Thomas Jefferson National Accelerator Facility (JLab), clustering aspects are studied 
in structure functions of deep inelastic lepton-nucleus scattering by using momentum distributions 
calculated in antisymmetrized (or fermionic) molecular dynamics (AMD) and also in a simple shell 
model for comparison. According to the AMD, the 9 Be nucleus consists of two a-like clusters 
with a surrounding neutron. The clustering produces high-momentum components in nuclear wave 
functions, which affects nuclear modifications of the structure functions. We investigated whether 
clustering features could appear in the structure function F2 of 9 Be along with studies for other light 
nuclei. We found that nuclear modifications of F% are similar in both AMD and shell models within 
our simple convolution description although there are slight differences in 9 Be. It indicates that the 
anomalous 9 Be result should be explained by a different mechanism from the nuclear binding and 
Fermi motion. If nuclear-modification slopes d(F^ / F®)/dx are shown by the maximum local den- 
sities, the 9 Be anomaly can be explained by the AMD picture, namely by the clustering structure, 
whereas it certainly cannot be described in the simple shell model. This fact suggests that the large 
nuclear modification in 9 Be should be explained by large densities in the clusters. For example, 
internal nucleon structure could be modified in the high-density clusters. The clustering aspect of 
nuclear structure functions is an unexplored topic which is interesting for future investigations. 

PACS numbers: 13.60.Hb, 13.60.-r, 24.85.+p, 25.30.-c 



I. INTRODUCTION 

Nuclear modifications of structure functions F% were 
found by the European Muon Collaboration (EMC) 
so that the phenomena is often called the EMC effect. 
Such modifications are now measured from relatively 
small x (~ 1CP 3 ) to large x (x ~ 0.8), where x is the 
Bjorken scaling variable. By using the data on nuclear 
structure functions, optimum parton distribution func- 
tions (PDFs) are proposed for nuclei [2j,l3j. Physics mech- 
anisms are different depending on the x region for pro- 
ducing the nuclear modifications. At small x, suppression 
of F2 occurs and it is known as nuclear shadowing. It is 
due to multiple scattering of a qq pair coming from the 
virtual photon. At medium and large x, modifications are 
understood by conventional models mainly with nuclear 
binding and Fermi motion of nucleons. However, it may 
not be possible to explain full experimental modifications 
by such mechanisms, which indicates that internal struc- 
ture of the nucleon could be also modified in a nuclear 
medium. For explanations of these physics mechanisms, 
the reader may look at Ref. [3]. 

In future, much details of the nuclear modifications 
will be investigated in lepton-nucleus deep inelastic scat- 
tering and hadron-hadron reactions. For example, a nu- 
clear modification difference between up and down va- 
lence quark distributions (u v and d v ) will be investigated 
by measuring cross sections of semi-inclusive ^ produc- 
tions [j| . It could lead to a possible solution of the long- 



standing NuTeV weak- mixing angle (sin Ow) anomaly Q 
from a viewpoint of the nuclear modification difference 
between u v and d v .7]. In addition, nuclear shell struc- 
ture of the EMC effect, so called "local EMC effect" @, 
will be investigated by measuring semi-inclusive reactions 
Q . There is also an issue of the nuclear- modification dif- 
ference between the structure functions of charged-lepton 
and neutrino reactions E| ■ It needs to be solved for 
a precise determination of nucleonic and nuclear PDFs. 
Such nuclear effects will be investigated by neutrino re- 
actions of the MINERVA project pj]. There will be also 
measurements at hadron facilities at RHIC (Relativistic 
Heavy Ion Collider), Fermilab (E906 experiment) [l3| . 
LHC (Large Hadron Collider), and possibly at J- PARC 
(Japan Proton Accelerator Research Complex) [14] . 

Measurements on the EMC effect at the Thomas Jef- 
ferson National Accelerator Facility (JLab) obtained an 
anomalous result for the beryllium-9 nucleus in compar- 
ison with measurements for other light nuclei 15]. It is 
anomalous in the sense that the magnitude of the nuclear- 
modification slope I d(F 2 A I ' F2) I 'dx | is much larger in 9 Be 
than the ones expected from its average nuclear density. 
From the experimental 9 Be radius, namely the average 
nuclear density, modifications of F2 at medium x are ex- 
pected to be much smaller than the ones of 12 C, whereas 
measured values are similar in magnitude. 

It is known in nuclear structure studies that the 9 Be 
nucleus has a typical cluster-like structure rather than 
a shell-like one 111, H/fll. It is like a cluster of two a 



( 4 He nucleus) particles with surrounding neutron clouds 
according to the studies of antisymmetrized molecular 
dynamics (AMD). This fact indicates that there exist 
higher- density regions than the ones expected from the 
average density by the shell model or the one estimated 
by the experimental charge radius. The high-density re- 
gions could contribute to larger nuclear modifications of 
the structure function F2. It could be a reason for the 
anomalous modification for 9 Be. Such a cluster structure 
could produce high-momentum components in the mo- 
mentum distribution of the nucleon, which is eventually 
reflected in modifications of quark momentum distribu- 
tions, namely the structure functions of nuclei. 

These considerations motivated us to investigate clus- 
ter aspects in the structure functions F^ for light nuclei, 
especially 9 Be. At this stage, there is no theoretical work 
on the nuclear-clustering aspect in high-energy nuclear 
processes, for example, in structure functions, although 
there are some studies on multi-quark clusters such as 
a six-quark state in 1980's. Our current studies are to- 
tally different from these works on multi-quark effects. 
In this article, we investigate possible nuclear clustering 
effects on the structure functions F^ within a convolu- 
tion model for describing F^ by using the AMD and 
shell-model wave functions. 

In. Sec. mi our theoretical formalism is provided for 
describing nuclear structure functions F^ ■ First, the con- 
volution model is introduced. Then, the AMD descrip- 
tion and a simple shell model are explained for calculat- 
ing nuclear wave functions. In Sec. IIII1 calculated nuclear 
densities are shown for He and 9 Be in order to illustrate 
the clustering structure in the 9 Be nucleus. Then, mo- 
mentum distributions are shown for these nuclei. The ra- 
tios F^/Ff are calculated in both AMD and shell mod- 
els, and they are compared with experimental data. Fi- 
nally, the nuclear modification slopes d(F^/F^)/dx are 
discussed. Our results are summarized in Sec. IIV1 



II. FORMALISM 

We explain a basic formalism for calculating the nu- 
clear structure functions F^ in the convolution approach 
together with antisymmetrized molecular dynamics and a 
simple shell model for calculating nuclear wave functions. 
These models are somewhat obvious within each commu- 
nity in structure-function and nuclear-cluster physicists. 
However, the following introductory explanations are in- 
tended that different communities could understand with 
each other. 



A. Nuclear structure functions in convolution 
approach 

The cross section of deep inelastic charged-lepton- 
nucleon (or nucleus) scattering is expressed by a lepton 
tensor multiplied by a hadron tensor W^ v : da ~ 



L^W^ [IGiHIJ. The hadron tensor is defined by 

W ^M = ^ J d^e^< <p\[MO,M0)]\p>, (i) 

where q is the virtual photon four- momentum, p is the 
momentum of the nucleon (or nucleus), and is the 
hadronic electromagnetic current. The hadron tensor 
Wfj, v is expressed by the imaginary part of the forward 
virtual Compton amplitude T„ v as Wau = -^ m (7) 1!/ )/(27r) 
by the optical theorem. 

The convolution model has been discussed in various 
articles within binding models for calculating nuclear 
structure functions, so that the detailed formalism should 
be found, for example, in Refs. [1, [l9j]. It indicates that 
a nuclear structure function is given by an integral of the 
nucleonic one convoluted with a momentum distribution 
of a nucleon in a nucleus as illustrated in Fig. [T] It is 
written in the hadron-tensor form as 

W* (pa ,q) = j d 4 PN S( PN ) W? v { PN , q) , (2) 

where pn and pa are momenta for the nucleon and nu- 
cleus, respectively, and S(j>n) is the spectral function 
which is the energy-momentum distribution of nucleons 
in the nucleus. The structure functions are generally ex- 
pressed in terms of two variables Q 2 and x defined by 



where is the nucleon mass, v is the energy transfer 
v = q° in the rest frame of a target nucleus, and q 2 is 
given by q 2 — (q ) 2 — q 2 ■ In the convolution picture 
of Eq. @, the process is described by two steps as il- 
lustrated in Fig. [U First, a nucleon is distributed in a 
nucleus according to the spectral function S(pn) with the 
nucleon momentum p^ , and then a quark is distributed 
with the momentum fraction x in the nucleon. The over- 
all quark momentum distribution is given by the convo- 
lution integral of these two distributions. 




FIG. 1: Convolution approach for nuclear structure functions. 
The 7*, q, N, and A indicate the virtual 7, quark, nucleon, 
and nucleus, respectively. A quark momentum distribution is 
described by the integral of a corresponding quark distribu- 
tion convoluted with a nucleon momentum distribution. 
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The hadron tensor W^ v is expressed in terms of two 
structure functions and W 2 as 



-W^{p A , q Yg^+Wi{p A ,q)^^, (4) 

Pa 



where g^ and p^ are defined by g^ u = g^ - q^q^/q 2 
and p^ = p^ — (p ■ q) q^/q 2 so as to satisfy the current 
conservation. The structure function F 2 is related to 
W 2 by F% = W a pa ■ q/MA, and its projection operator 



is given by |2l|, |2j 



r 2 



M A pa • q 



2pi 



Pi 



(5) 



which satisfies P^W^, 



Fo- The mass of the nucleus 



is denoted by M4. Applying the proj ection operator on 
both side of Eq. ©, we obtain jJHJ ^ 



F 2 A (x,Q 2 ) 



dyf(y)F 2 N (x/y,Q 2 ), 



(6) 



where i 7 ^ 4 and i 7 ^ are structure functions for the nucleus 
and nucleon, and y is the momentum fraction 



y 



M A p N -q Ap-jy 



M N p A ■ q 



Pa 



(7) 



where p + is a light-cone momentum [p + = (j>° +p 3 )/\/2]. 
It should be noted that the upper bound of the variables 
x and y is A for nuclei. The function f(y) indicates a 
light-cone momentum distribution for the nucleon, and 
it is given by 



d 3 PN y8[y 



pn ■ q 



where rii is the number of the nucleon in the quantum 
state i, and the summation is taken over the occupied 
states. Here, the spectral function is given by 



S(Pn) = -J^2 n i \<l>i(j>N)\ 2 S (p% 



M A 



M A 2 _ 



Pn 



(9) 

where M A -i is the mass of residual one-hole state, and 
4>i(PN) is the wave function of the nucleon. Here, p^/Q 2 - 
type higher-twist effects [2l[ are not included in the con- 
volution equation. The function f(y) is normalized so as 
to satisfy the baryon-number conservation J Q dyf(y) = 1 
by taking / dp%yS{p N ) = 1 @, [!|. 

The wave functions of the nucleon are calculated non- 
relativistically, and then they are used for the relativistic 
description in obtaining light-cone distributions by Eq. 
©. It could lead to an issue of normalizing the non- 
relativistic wave function because there is no solid rela- 
tivistic framework to use the non-relativistic functions. 
Here, the wave functions are normalized to satisfy the 
condition J dyf(y) — 1, where there is an extra factor 



of p° n /Mm in front of |</>(pjv)| 2 . As noticed in the third 
article of Rcf . [25[ , this factor does not appear if a mass 
factor (Mn/p%) is included in the convolution formal- 
ism. However, such an overall normalization difference 
does not affect our results in Sec. Mil 
The separation energy Ei is defined by 



e t = (M A -i + M N ) - Ma- 



(10) 



It is the energy required to remove a nucleon from the 
state i. In our actual calculation, we average over all 
the nucleons for estimating the average separation en- 
ergy (sj— >< £ >). If a non-relativistic approximation is 

applied for the expression jM~ A \ i +p^, p% and < e > 

are related by considering the S function for the energy 
conservation as 



p% = M N - < e > - 



Pn 



2M A -i ' 



(11) 



where Ma-% is replaced by Ma-i (the ground-state mass 
of the A— 1 nucleus). It should be noted that the residual 
nucleus A — i could be in an excited state and that many- 
body breakup processes could be also possible in the final 
state. Therefore, the separation energy is, in general, not 
a simple difference between the two nuclear binding ener- 
gies in the initial and final states, since the final nucleus 
would not be in the ground state. It means that theoret- 
ical separation energies depend how they are estimated. 
For example, they vary depen ding whether models in- 
clude short-range correlations 24, [26[ and many-body 
breakup processes (27j . In our work, experimental sepa- 
ration energies are taken from (e, e'p) and (p, 2p) exper- 
iments. 

Equation ^ indicates that the nuclear structure func- 
tion F 2 is split into two parts: the light-cone momen- 
tum distribution of the nucleon and the nucleonic struc- 
ture function F 2 . If there is no nuclear medium effect 
on the nucleonic structure function F 2 N (x,Q 2 ), nuclear 
modifications should come solely from the nucleonic dis- 
tribution part, which contains the information on nuclear 
binding and Fermi motion of nucleons. These effects are 
reflected in the light-cone momentum distribution of Eq. 
(|8]), namely in the momentum distribution of the nucleon 
and the energy-conserving S function. For calculating the 
distribution f(y), we need a realistic model for the wave 
function c/)(pn). In our work, we calculate it in two theo- 
retical models: an antisymmetrized molecular dynamics 
and a simple shell model. They are introduced in Sees. 
iHBlandlrrCl 



B. Antisymmetrized molecular dynamics 

This work is intended to investigate a possible clus- 
tering effect on the structure functions of deep inelastic 
scattering (DIS). There is a theoretical method, antisym- 
metrized molecular dynamics (AMD) [28| or fermionic 
molecular dynamics (FMD) [29J, which is developed for 
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describing clustering aspects of nuclei as well as shell- 
like structure on an equal footing. Hereafter, we use the 
nomenclature AMD for this theoretical method. 

There are nuclei which exhibit density distributions 
of separate clusters. For example, the 8 Be nucleus has 
two separate peaks, which correspond to two a nuclei, 
in its density distribution according to a Monte Carlo 
calculation for the eight-body system by using realistic 
NN (N: nucleon) potentials [3Qj |. It suggests that some 
nuclei tend to form a-like clusters within their structure 
since the a is a tightly bound nucleus. 

A simple and yet very useful and consistent theoretical 
method is provided by the AMD method. The AMD has 
a number of advantages, for example, that there is no 
assumption on nuclear structure, namely shell or cluster 
like configuration, and that simple and systematic studies 
are possible from light to medium-size nuclei. A nuclear 
wave function is given by the Slater determinant of single- 
particle wave packets: 



1 



det[v5i(r*i), ^2(^*2), • • ^ipa^a)] 

(12) 



VAl 

¥>2(ri) ip 2 (r 2 ) 



^2{ta) 



Wa{ta) 



Here, a nucleon is described by the single-particle wave 
function 



(13) 



where \i an d t, indicate spin and isospin states, respec- 
tively. The function <f>i(rj) is the space part of the wave 
function, and it is assumed to be given by the Gaussian 
functional form: 



3/4 



exp 



(14) 



where v is a parameter to express the extent of the wave 
packet. The center of the wave packet is given by Zij 'y/u. 
We should note that Zi is a complex variational param- 
eter. Its real and imaginary parts indicate nucleon posi- 
tion and momentum, respectively [28j : 



< r\(f>j > _ ReZj < (pi\p\<fH > 

< 4>i\4>i > y/V ' < (j)i\4>i > 



2y/UlmZ t . (15) 



A nuclear state is an eigenstate of the parity, so that the 
following parity-projected wave function is used: 

|4 ,± (ri,r 2 ,- • ^Ta)) =^=[ |$(ri,r 2 ,- • ;f A )) 

±\$(-r 1 ,-r 2 ,--;-r A )) ]■ (16) 



As for the NN interactions, we use the following po- 
tentials: 

2- body: V 2 = (1 - m - mP a P T ) 

x v 21 e- {r/r2l)2 +v 22 e~ {r/r22) 

3- body: V 3 = v 3 5 3 (f 1 - r 2 ) <5 3 (r 2 - r 3 ), 

LS: V LS = v LS 



g-(r/r L si) _ e -(r/r LS 2Y 



x P( 3 0)L-S, 



(17) 



where m, v 21 , v 22 , r 21 , r 22 , v 3 , v LS , r LS1 , and r LS2 are 
constants. The two-body interaction part mP„P T indi- 
cates the Majorana term with spin and isospin exchange 
operators (P a , P T ). The three-body part is a contact in- 
teraction form, and P( 3 0) is the projection operator of 
the triplet-odd ( 3 0) state (spin 5=1, ang ular momentum 
L=odd) in the two- nucleon system [31(. The Coulomb 
interaction is also considered in our analysis. The con- 
stants m, V21, • • • are taken from Ref. [1 €31 ] except for 
v 3 and VlSi which are fixed so as to reproduce binding 
energies of considered nuclei under the radius constraint: 
v LS = 2000 MeV, v 3 = 4000, 3300, 2000 MeV for 4 He, 
9 Be, 12 C, respectively. Here, we should be careful to take 
into account the effect of center-of-mass motion [32| . 

The AMD wave functions contain the parameters Zi 
and v, which are determined by minimizing the system 
energy with a frictional-cooling method. Time develop- 
ment of Zi is described by the time-dependent variational 
principle: 



12 



dt 



;$(z) 



' dt 



H\HZ)\ 



'MZ)\${zy, 



(18) 



It leads to the equation of motion. Introducing two ar- 
bitrary parameters A and \x for practically solving the 
equation of motion, we obtain 



l -Zi = {X + l ^) — 



(19) 



Here, \i is a friction parameter which should be a negative 
number. By solving this equation, the parameters Zi are 
obtained. From the obtained parameters, the densities 
in coordinate and momentum spaces are calculated by 



p(r) = 



Pip) = 



— V 

2-kv 1 ^ 



exp 



(20) 



(21) 



where Bij = J dfip\(r)ipj(r). This momentum distri- 
bution is used for calculating the light-cone momentum 
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FIG. 2: (Color online) Shell levels for 3 He, 4 He, 9 Be, and 12 C. 

distribution in Eq. Then, using the convolution 

equation of ©, we obtain the nuclear structure func- 
tions, which include clustering effects described by the 
AMD. 



take the constant u> for each nucleus so as to obtain the 
same matter radius calculated by the AMD. 

D. Deuteron wave function 

Experimental data are listed by ratios F^ixj/F®, 
where is the structure function of the deuteron, for 
showing nuclear modifications in the structure functions. 
Since the deuteron is a bound two-nucleon system, a com- 
mon wave function is used in Eq. © for calculating the 
structure function F^(x,Q 2 ) of Eq. © in both AMD- 
and shell-model analyses. Here, we take the deuteron 
wave function given by the Bonn group in Ref. (36j . 



C. Simple shell model 

In order to compare with the AMD results at this 
stage, we also calculate the nuclear spectral function by 
using a simple shell model, because the current wave 
functions and NN interactions in AMD are simple Gaus- 
sian forms. If much detailed studies become necessary in 
future, we may consider to use more sophisticated mod- 
els, for example, a density-dependent Hartree-Fock [33j 
or a detailed shell model such as NuShell (OXBASH) 
34 1. 

As a shell model, we take a simple harmonic oscilla- 
tor model. Nucleons are assumed to move in an aver- 
age central potential created by interactions of all the 
nucleons in a nucleus. Then, the nucleons are treated 
independently with each other. A simple and yet realis- 
tic choice of the potential is the harmonic-oscillator type 
(M/vw 2 r 2 /2). Its wave function is separated into radial- 
and angular-dependent parts: 



(22) 



where r, 9, and <f> are spherical coordinates, and n, £, and 
m are radial, azimuthal, and magnetic quantum num- 
bers, respectively. The function Yi m (9, <f>) is the spherical 
harmonics, and the radial wave function is given by [3o| 



Rni(r) 



1 2n 2i + 3 (n-iy. 
\T(n + e+1/2)] 3 



r e 



where L^_}( 2 (x) is the Laguerre polynomial, and K is 
defined by K = \/Mnuj. 

In the following analysis, the light nuclei, 3 He, 4 He, 
9 Be, and 12 C are considered, so that the low-energy lev- 
els, lsi/2 and lps/2, are taken into account as shown in 
Fig. [2] The only parameter in the model is w, which 
is fixed by a nuclear radius. The constants of the AMD 
model are determined so as to explain experimental nu- 
clear charge radii. Then, nuclear matter radii are calcu- 
lated by using obtained AMD densities. Since it is the 
purpose of this work to investigate a difference between 
the structure functions of the AMD and shell models, we 



(23) 



E. Experimental information on separation 
energies, binding energies, and charge radii 

In calculating the structure functions, experimental in- 
formation is needed for separation energies, binding en- 
ergies, and charge radii. The binding energies are taken 
from Ref. [37j and they are listed in Table HI They are 
used for calculating nuclear mass: Ma = Z M p + N M n — 
B, where M p and M n are proton and neutron masses, Z 
and N are atomic and neutron numbers, and B is the 
binding energy. 

Experimental nuclear charge r.m.s. (root-mean- 
square) radii are listed for the deuteron (38|, 3 He [39j j . 
4 He H3], 9 Be [U, and 12 C [42] in Table H Using these 
charge radii and binding energies, the constants in the 
AMD model are determined. The matter r.m.s. radii are 
then calculated in the AMD by using Eq. (l20l) . There 
are slight differences between V '< r' 1 > c and \/< r 2 > m 
in the AMD for 9 Be and 12 C. This is due to the effect of 
Coulomb force. 

The separation ener gies are taken from experimental 
measurements for 4 He [43l| . 9 Be and 12 C and they 
are listed in Table fl] A theoretical estimate 11.4 MeV is 
listed just for information because there is no available 
data for 3 He. It was obtained by using a spectral function 
calculated by the Faddeev method with the Reid soft-core 
potential [45j j . 



TABLE I: Experimental data for mean separation energies, 
binding energies per nucleon, and charge root-mean-square 
radii. The matter radii are calculated in the AMD model 
except for the deuteron so that charge radii agree with the 
data. The asterisk * indicates a theoretical estimate of Ref. 
[45( ] because the experimental data is not available. 



Nucleus 


< £ > 

(MeV) 


B/A 
(MeV) 


V< r 2 > c 
(fm) 


V< r 2 > m 
(fm) 


D 


2.22 


1.11 


2.10 


2.10 


3 He 


(11.4*) 


2.57 


1.96 


1.96 


4 He 


20.4 


7.07 


1.68 


1.68 


9 Be 


24.4 


6.46 


2.52 


2.61 


12 C 


22.6 


7.68 


2.47 


2.48 



5 



It should be also noted [46| that the separation energy 
20.4 MeV of 4 He [43| was obtained by using the data only 
in the peak region of the energy spectrum of 4 He(p,2p) 3 H 
and a continuum region is not included. The separa- 
tion energy should be calculated by the average energy 
weighted by the spectral function: 

< £ >= J dE N d 3 p N E N S(E N ,p N ). (24) 

We notice that theoretical estimates are usually larger 
than this value (20.4 MeV) for 4 He [H,[27j]. For exam- 
ple, 28.2 MeV is obtained in Ref. [26j |. where the average 
kinetic energy estimated by the ATMS (Amalgamation of 
Two-body correlation functions into Multiple Scattering 
process) method is employed, and then the Koltun sum 
rule is used for estimating the separation energy < e >: 
B/A= [<e> - <T> (A-2)/(A- l)]/2, where B/A 
is the binding energy per nucleon and < T > is the 
average kinetic energy. However, it is very difficult to 
calculate a reliable value of the separation energy. The 
experimental separation energies < s > were obtained in 
nucleon-knockout reactions by observing peaks of single- 
particle excitations and they do not include the contri- 
bution from continuum states of the residual nucleus. 
Therefore, the mean separation energies would be un- 
derestimated. In this work, we estimated the clustering 
effect without the continuum, which needs to be consid- 
ered in future for detailed comparison with data. 



the polar and azimuthal angles 9 and <f>. Although the 
4 He densities are same in both models, they are differ- 
ent in 9 Be. Since the angular integrals have been done, 
the cluster structure is no longer apparent in the AMD 
density of 9 Be in Fig. SJ However, the cluster effects 
are reflected in the slightly larger densities at r ^2 fm 
and the depletion at r = 0, due to the existence of two 
separate clusters. 

Instead of the coordinate-space density, the 
momentum-space density |0(p*/v)| 2 is used for cal- 
culating the light-cone momentum distribution by Eq. 
((SJ). Calculated momentum-space densities are shown in 
Fig. [5] for the nuclei 4 He and 9 Be. We explained in Sec. 
Ill El that the same radii are taken in both AMD and 
shell models. As a result, both momentum distributions 
of 4 He are almost the same. However, the distributions 
are much different in 9 Be. It is important to find that the 
momentum distribution of the AMD is shifted toward the 
high-momentum region in 9 Be because of the clustering 
structure. This is caused by the fact that the dense 
regions, namely the two clusters, are formed within the 
9 Be nucleus. If nucleons are confined in the small space 
regions of the clusters, it leads to an increase of high 
momentum components, which is clearly shown in Fig. 

m 

Now, using the obtained momentum distributions to- 
gether with Eqs. ©, ©, and ©, we calculate the nu- 
clear structure functions. The structure function of 9 Be is 



III. RESULTS 

First, nuclear densities are shown in the AMD model. 
The focused nucleus is 9 Be for investigating the anoma- 
lous EMC effect in the structure function F2, so that 
coordinate-space densities are shown in Fig. [3] for this 
9 Be nucleus as well as 4 He, as an example, for compar- 
ison. It is interesting to find two density peaks within 
the 9 Be nucleus, whereas the 4 He density is a monotonic 
distribution. As mentioned in Sec. UH it is the advantage 
of the AMD method that it docs not assume any specific 
structure, cluster- or shell-like configuration, on nuclei. 
The 4 He is a tightly bound nucleus and it is well described 
by the usual shell-like structure, which is judged by the 
monotonic density distribution in Fig. [31 However, the 
situation is apparently different in 9 Be. The figure sug- 
gests that two dense regions exist in 9 Be although such 
a phenomenon does not exist in the shell model. It indi- 
cates that the 9 Be nucleus consists of two a-like clusters 
with surrounding neutron clouds. This clustering could 
produce different nuclear medium effects from the ones 
expected by the shell model. In particular, it could in- 
fluence the nucleon momentum distribution, eventually 
quark momentum distributions, within the 9 Be nucleus. 
Furthermore, dense regions could alter the internal struc- 
ture of the nucleon. 

Next, coordinate-space densities are compared in both 
AMD and shell models in Fig. H]by taking averages over 




FIG. 3: (Color online) Upper and lower figures indicate 
coordinate-space densities of 4 He and 9 Be, respectively, cal- 
culated by the AMD. Here, the densities are shown by taking 
integrals over the coordinate z: f dzp(x,y, z). 
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FIG. 4: (Color online) Coordinate-space densities are shown 
for 9 Be ( 4 He) in the AMD and shell models by the wide-solid 
(narrow-solid) and dashed (dotted) curves, respectively. The 
densities are integrated over the angles 9 and (j> for showing 
the curves in this figure. The clustering structure in the AMD 
gives rise to a modification of the density distribution in 9 Be, 
whereas both densities are the same in 4 He. 
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FIG. 6: (Color online) Theoretical structure-function ratios 

Fi He /FP are compared with experimental data of SLAC- 
E139 NMC H, and JLab The solid and dashed 

curves indicate AMD and shell model results, respectively, 
calculated at Q 2 =5 GeV 2 ; however, both curves overlap each 
other. The experimental data are taken at various Q 2 points. 
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FIG. 5: (Color online) Momentum-space densities are shown 
for 4 He and 9 Be in the AMD and shell models. The clustering 
structure in 9 Be gives rise to an excess of high-momentum 
components in the AMD. 
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FIG. 7: (Color online) Comparison with JLab experimental 
data of F 2 B£ /F,p. Notations are the same as the ones in Fig. 
[6] The differences between the solid (AMD) and dashed (shell 
model) curves are now clear in 9 Be. Isoscalar corrections are 
applied as explained in the main text. 



shown together with the one of 4 He as an example of non- 
cluster-like nuclei in order to illustrate clustering effects 
on the structure function of 9 Be. In Figs. [6] and [71 our 
theoretical ratios F 2 He /F 2 D and F 2 Be /F® are compared 
with the available experimental data of the SLAC (Stan- 
ford Linear Accelerator Center)-E139 E^, NMC (New 
Muon Collaboration) [Hj], and JLab fisjj. The AMD 
and shell-model ratios are shown by the solid and dashed 
curves, respectively, and they are calculated at a fixed 
Q 2 point (Q 2 = 5 GeV 2 ). Experimental data are taken 
at various Q 2 points, and only the data with Q 2 > 1 
GeV 2 are shown in Figs. [6] and [Jj The JLab measure- 
ments include the data with small invariant mass W, 
where the process is not considered to be deep inelastic. 
Therefore, the data with W 2 < 3 GeV 2 are shown by 
the open circles. In showing the ratios of non-isoscalar 
(Z ^ N) nuclei, isoscalar corrections are applied in Ref. 
15] by including smearing corrections. Since we can- 
not access to the specific smearing corrections in the 



JLab analysis, we simply used the isoscalar corrections 
[(Ff + F 2 n )/2]/[(ZF| + NF£)/A], where if and F 2 ™ are 
the structure functions of the proton and neutron, re- 
spectively, by using the PDFs of the MSTW08 [111 in 
the leading-order (LO) of a s . We have checked that our 
corrections are almost the same as the corrections in the 
JLab analysis in 9 Be fl5| . 

From Figs. [6]and[7J we find that our theoretical ratios 
have a tendency consistent with the data in the sense 
that the ratio decreases at medium x and it increases at 
large x. These decrease and increase are caused by the 
nuclear binding and the nucleon's Fermi motion, respec- 
tively, in our convolution picture. However, it is also clear 
that the simple convolution description is not sufficient to 
explain the whole experimental nuclear modifications be- 
cause there are differences between the theoretical curves 
and the data. 

There are two major reasons for the differences. First, 
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short-range nucleon-nucleon correlations have not been 
included in calculating the spectral function [24|, [26j . 
They change the theoretical ratios toward the experimen- 
tal data at x = 0.6 — 0.8. The purpose of our studies is 
to investigate whether or not a possible clustering signa- 
ture appears in deep inelastic lepton-nucleus scattering. 
Since this is the first attempt to investigate the cluster 
effects, we did not include such an effect. In future, we 
may consider to study more details. 

Second, there could be a modification of nucleon itself 
inside a nuclear medium. As explained in Refs. [3. [l9j. 
such a nucleon modification was originally proposed as a 
Q 2 rescaling model. Nucleons could overlap in a nucleus 
since the average nucleon separation and nucleon diame- 
ter are almost the same. The overlap then gives rise to a 
confinement radius change for quarks, which appears as 
a modification of quark momentum distribution, namely 
a modification of the structure function F-i- A possible 
internal nucleon modification was investigated in Refs. 
[49l |50] ] in comparison with the data. Since it is not the 
purpose of this work to step into such details, especially 
in comparison with the data, we leave it for our possible 
future studies. 

The anomalous data was reported for 9 Be by the JLab 
experiment [l5[ by taking a slope of the ratio F^/F® 
with respect to the Bjorken variable x in the region 
0.35 < x < 0.7. As shown in Fig. 4 of Ref. [Hj], the mag- 
nitude of the 9 Be slope is too large to be expected from 
its average nuclear density in comparison with the ones 
of other light nuclei. We calculate corresponding the- 
oretical slopes by taking the derivatives d(F^ / F®)/dx 
(= cLRemc /dx) at x = 0.35. The JLab data are plotted 
by the average density calculated by a Greens Function 
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FIG. 8: (Color online) Comparison with JLab data on the 
slope \<IRemc /dx\, where Remc = F% /F} ' ■ The open circles 
and crosses are the theoretical slopes calculated by the AMD 
and shell models, respectively. The JLab data are shown 
by the filled circles with errors. The abscissa is the aver- 
age density p defined by A/[4:TY < r' 2 >^J 2 /3] with the charge 
r.m.s. radius %/<r 2 > c - The curve indicates a smooth func- 
tion I dREMc/dx I = a(p— p D ) b to fit the JLab experimental 
data except for 9 Be. 



Monte Carlo method [1TJ with the multiplication factor 
(A — 1) I A for removing the struck nucleon. 

Although this theoretical density estimate would be 
reliable, we first show the slope by a purely exper- 
imental quantity by defining an average density as 

A/[4ir(r 2 ) 3 c /2 /3], where VkW > c is the experimental 
charge r.m.s. radius, in Fig. [5] instead of the specific 
theoretical density. The experimental charge radii are 
taken from Table |TJ The theoretical slopes in the AMD 
and shell models are shown by the open circles (O) and 
crosses (x), respectively, for 4 He, 9 Be, and 12 C. Since 
the experimental separation energy is not available for 
3 He, the theoretical slopes are not calculated for 3 He. 
The JLab data are shown by the filled circles with er- 
rors. In order to illustrate how the 9 Be slope deviates 
from the other nuclear ones, a curve is given in Fig. 
[5] by fitting the data without the 9 Be data in a sim- 
ple functional form, | dRsMc/dx \ = a(p — p D ) b where 
p = A/ [Air < r 2 > Z J 2 /3] and p D is the density of the 
deuteron. The parameters a and b are determined and 
we obtain | dREMC I dx \ — 1.35(p — P D ) ' 906 , which is the 
curve in Fig. [H It is obvious that the 9 Be slope is anoma- 
lous in the sense that the data significantly deviates from 
the curve. 

The magnitudes of the theoretical slopes are rather 
small in comparison with the data, and they are about 
half or less of the experimental ones in Fig. [5] This was 
already obvious from Figs. [5] and [7] that the magnitudes 
of the theoretical slopes are smaller than the experimen- 
tal ones. As explained, the differences could be caused by 
the short-range correlations and internal nucleon modifi- 
cations. An interesting result is that the clustering effects 
are not apparent in the slope by looking at both AMD 
and shell-model results for 9 Be although there are some 
differences in the momentum distributions of Fig. [5] and 
in the structure functions of Fig. [7] In other nuclei, 
both theoretical slopes are almost identical. This is un- 
derstood in the following way. In the medium- a; region, 
the nuclear modifications can be described mainly by the 
first two moments of the nucleon momentum distribution 
f(y). These moments are expressed by the av erag e sepa- 
ration and kinetic energies, < e > and < T > [25], which 
are similar in both models. It leads to the small differ- 
ences between the AMD and shell models in the slope 
dR E Mc/dx. 

The small difference between the AMD and shell-model 
slopes in 9 Be suggests us to look for another reason to 
explain the anomalous JLab data. As we noticed in Fig. 
131 the high density regions are created locally in 9 Be ac- 
cording to the AMD model. The higher densities could 
contribute to extra nuclear modifications in the struc- 
ture function F^ by an additional mechanism which is 
not considered in our simple convolution picture. In or- 
der to find such a possibility, we plot the same slope by 
taking the maximum local density as the abscissa. The 
maximum density p m ax, of course, depends on theoretical 
models to describe the nuclei. The maximum positions 
are located at r — for 4 He in both AMD and shell mod- 
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els and also for 9 Be in the shell model. However, they 
are at different points (r ^ 0) for 9 Be in the AMD due 
to the cluster structure and for 12 C in both models due 
to a p-wave contribution to the density. 

The slopes are shown in Fig. |9]by taking the maximum 
local density p m ax as the abscissa. The maximum densi- 
ties are almost the same in 3 He and 4 He, so that they are 
plotted at the same position of p max . However, they dif- 
fer in 9 Be and 12 C. Although the difference between p max 
(AMD) and p m ax (Shell) in 12 C is not as large as the one 
in 9 Be, it seems that there exist some clustering effects 
also in 12 C. The JLab data of 9 Be and 12 C and the the- 
oretical slopes are plotted at two different density points 
of the AMD and shell models. The curve indicates a fit 
to the JLab data with the 12 C data at the shell- model 
density point by excluding the 9 Be data. It is given by 
I dRsMc/dx | = 0.82l(p max - p m ax D ) ' 646 , where p maXD 
is the maximum density for the deuteron. 

In the usual convolution calculation, the spectral func- 
tion is given by the averaged nuclear density distribu- 
tion, and thus the inhomogeneity of the nuclear density 
is washed out. In such calculation, the average nuclear 
density of 9 Be is lower than that of 12 C or 4 He as shown in 
Fig. [51 which is the origin of the "anomalous" EMC ratio 
of 9 Be observed at the JLab. However, the maximum lo- 
cal density of 9 Be is, as shown in Fig. higher than that 
the ones of 4 He and 12 C, and the EMC ratio of 9 Be can be 
treated "normally" . Here, the "normal" means that the 
result of 9 Be is consistent with the smooth curve deter- 
mined by the EMC results of other nuclei, and thus 9 Be 
does not have anomalous dependence on the nuclear den- 
sity anymore. In this sense, the anomalous 9 Be result is 
"explained" as a normal one by the maximum local den- 
sity; however, it does not mean that physics mechanism 
is clarified. 
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FIG. 9: (Color online) Comparison with JLab data on 
\cLRemc / dx\. The abscissa is the maximum local densities 
theoretically calculated in the AMD and shell models. Two 
models produce different densities in 9 Be and 12 C, so that 
the JLab and theoretical slopes are plotted at different den- 
sity points. The curve indicates a smooth function to fit the 
JLab experimental data except for 9 Be with the 12 C data at 
the shell-model density. 



What we emphasize here is that the EMC ratio or the 
nuclear structure function itself could consist of the mean 
conventional part and the remaining one depending on 
the maximum local density. The remaining part is surely 
associated with the inhomogeneity of the nuclear density, 
before taking the average of nuclear wave function, given 
by the nuclear cluster structure. It could be nuclear- 
medium modification of the nucleonic structure function. 
It is well known that the cluster structure is well devel- 
oped in Be and that the light mass region with A < 20 is 
very suitable to study the cluster structure. Therefore, 
it is reasonable that, although the cluster-structure effect 
in the EMC ratio is not seen in the medium and large nu- 
clei so far, we can now observe it in the beryllium isotope 
region. Such cluster structure will be investigated in the 
light- mass region by future JLab experiments [521 ] - 

It is interesting to find that the "anomalous" JLab 
data for 9 Be can be explained if it is plotted by the max- 
imum local density at the cluster positions because the 
9 Be (AMD) data is very close to the curve in Fig. [9] 
On the other hand, the 9 Be data remains anonymous if 
it is plotted by the shell-model density because the 9 Be 
(Shell) data significantly deviates from the curve. Such a 
tendency also exists in 12 C but it is not as serious as the 
9 Be case. If the average nuclear density of 9 Be is used 
in showing the slope data, the clustering effects are not 
clearly reflected. Here, it is important to point out that 
the 9 Be data agrees with the other nuclear data if they 
are plotted as a function of the maximum density. This 
fact implies that the physics mechanism associated with 
the high densities, for example due to the clusters in 9 Be, 
could be the origin for explaining the nuclear-modification 
slopes of the structure functions. One of the possi- 
ble mechanisms is the modification of internal nucleon 
structure caused by nuclear medium effects at the high 
density regions. 



IV. SUMMARY 

Nuclear modifications of structure function Fi were in- 
vestigated for finding a possible signature of clustering 
structure in nuclei. The convolution model was used for 
describing nuclear structure functions, where momentum 
distributions of the nucleon were calculated in the AMD 
and shell models. According to the AMD, the 9 Be nu- 
cleus has a clear clustering structure of two a-like clus- 
ters with a surrounding neutron. Because of the cluster 
formation in 9 Be, high-momentum components increase 
in the nuclear wave function of the AMD in comparison 
with the distribution of the shell-model one. Although 
there are some differences between the structure func- 
tions of 9 Be in the AMD and shell models, the differ- 
ences are rather small in our simple convolution picture. 
Therefore, an anomalous EMC effect found for 9 Be at 
JLab should come from other effects such as the inter- 
nal nucleon modification due to the high-density regions 
created by the clustering. 
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The following points are the major results in this work: 

(1) For the first time, the nuclear structure functions 
F^ are calculated in a model with clustering struc- 
ture in nuclei. Then, they are compared with the 
structure functions of the shell model to clarify the 
clustering effects. 

(2) The clustering configuration in the 9 Be nucleus pro- 
duces high-momentum components in the nuclear 
wave function. It leads to a modification of the 
light-cone momentum distribution for nucleons in 
the 9 Be nucleus. 

(3) Because of the high-momentum components due to 
the cluster formation, the nuclear structure func- 
tions F^ are modified; however, the modifications 
are not very large within the simple convolution 
description. 

(4) The anomalously large nuclear effect for the slope 
| d(F^ / F^)/dx | of the 9 Be nucleus observed at 
JLab can be explained if the slope is plotted by 
the maximum local density calculated in the theo- 
retical model (AMD) with clustering structure. 



(5) Since the nuclear-modification slopes are explained 
by the maximum densities of nuclei, the physics 
mechanism of the anomalous nuclear effect could 
be associated with the high densities in the clus- 
ters of 9 Be. This fact implies that internal nucleon 
modifications due to the high densities could be 
the origin of the 9 Be anomaly, although careful es- 
timations should be made on effects of short-range 
nucleon-nucleon correlations. 

This work is the first attempt to connect the DIS struc- 
ture functions to the clustering structure in nuclei. This 
kind of research field is an unexplored area, and further 
theoretical studies are needed for clarifying clustering ef- 
fects in the structure functions. 
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